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ABSTRACT 

We examine the evolution of the zonal flow pattern in the upper solar convection 
zone during the current extended solar minimum, and compare it with that during 
the previous minimum. The results suggest that a configuration matching that at the 
previous minimum was reached during 2008, but that the flow band corresponding to 
the new cycle has been moving more slowly towards the equator than was observed in 
the previous cycle, resulting in a gradual increase in the apparent length of the cycle 
during the 2007 - 2008 period. The current position of the lower-latitude fast-rotating 
belt corresponds to that seen around the onset of activity in the previous cycle. 

Subject headings: Sun: activity, Sun: helioseismology, Sun: rotation 

1. INTRODUCTION 

The onset of Solar Cycle 24 appears to be later than expected; nearly thirteen years after the 
Cycle 23 minimum in mid- 1996 there is still very little surface magnetic activity. In this Letter we 
examine the large-scale zonal flow pattern seen in the upper convection zone in recent observations 
and compare the results with those for the previous solar minimum. 

The migrating zonal flow pattern known as the torsional oscillation was first detected by 



Howard & Labontel (119801) in Doppler measurements at the solar surface carried out at the Mount 
Wilson Observatory; it consists of belts of slightly faster than average rotation that migrate from 
mid-latitudes to the equator and poles. The migration of the zonal flow bands during the solar 
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cycle is closely connected to the migration of the magneti c activity belt. The flow pattern was first 
detected helioseismically by iKosovichev & Schoul (|1997|) in early /-mode data from the Michel- 
son Doppler Imag er (MDI) abo ard the Solar and Heliospheric Observatory spacecraft, and was 
seen to migrate by lSchoul (|1999|) . Using p-mode data it is possible to resolve the depth penetration 
of the flows; however, only the North -South symmetric part of the pattern can be detected using 
global-mode helioseismology. The first four/five years of zonal-flow results from G lobal Oscil 



lation Networ k Group (GONG) and M DI p-mode and /-mode data were reported by iHowe et al. 



(120001) . while lAntia & Basul (120011) and ISchoul (12003b drew attention to the high -latitude part of 



the pa tt ern. The further evo l ution of the patter n in Cy cle 23 was describ ed by IVorontsov et al 



(120021) : lBasu& Antial (120031) : IHowe et all (120051) . and by IHowe et all (120061). who also co nsidered 
the updated Doppler observations from Mount Wilson. The analysis of lAntia et al.l (|2008) extends 
these investigations into the current solar minimum. It became evident that t he torsional o scilla - 
tion pattern penetrates through a significant fraction of the convection zone; IHowe et al.1 (|2005l) 
also found some phase variation with depth in the lower-latitude, equatorward-moving part of the 
pattern, in the sense that the increase in rotation rate happens at earlier times at deeper depths. The 
lower-latitude branch of the pattern migrates equatorward with the activity belts, starting at about 
45 deg latitude a few years before the previous solar minimum and reaching the equator a couple 
of years after maximum; the broader, stronger belt poleward of 45 deg propagates poleward on 
a rather shorter timescale. The whole phenomenon is usually believed to be a side-effect o f the 
dynamo that drives the solar cycle, either directly throug h the Lorent z effect (|Schuesslerll 198 1|) or 
indirectly through geostrophic effects in the activity belt dSpruiJ2003l) . Nevertheless, the bands are 
clearly detectable even when there are few or no organized active regions, as was seen in the early 
years of Cycle 23 and again at the dawn of Cycle 24. The magnetic or sunspot butterfly diagram 
does not exhibit purely periodic behavior; similarly, there is no reason to expect that the zonal flow 
pattern should repeat precisely from cycle to cycle, but it is interesting to compare the flows seen 
at different minima and the way they evolve. 

In this Letter, we consider the behavior of the flows around the current and previous solar 
minima. We concentrate on the equatorward-moving part, basing our comparisons on latitudes 
equatorwards of 45 deg, and on the upper part of the convection zone where the signal is clearest. 



2. DATA 

We now have inferences of the solar interior rotation rate, Q(r, 6) where r is radius and 
9 is latitude, covering the period from mid- 1995 (GONG) and mid- 1996 (MDI) to early 2009. 
These inferences were obtained using 2-dimensional regularized least-squares inversions on sets 
of medium-degree p-mode data from a total of 137 GONG time series (108 days, with start times 
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at 36-day intervals) and 64 non-overlapping MDI time series of 72 days each, contiguous except 
for the interruptions due to SQH O problems in 1998 and 1999. Details of the data analysis and 
inversion techniques are given in iHowe et al.1 ([2005) and references therein. The inversions are 
carried out on a mesh with 48 equal intervals in latitude and 50 non-uniform intervals in radius; 
the radial mesh was chosen to be approximately uniform in acoustical radius. To reveal the tor- 
sional oscillation signal, at each location a temporal mean over all epochs was subtracted from 
the rotation-rate inference. The mean was calculated separately for MDI and GONG inversions, 
before the residuals were combined into one time series. 

The synoptic magnetic index used here has been derived from a combination of KPVlfl and 
SOUS VSMlmagnetograms. 



3. RESULTS AND ANALYSIS 

Figure \T\ shows the zonal flow pattern as a function of latitude and time at O.99_R . The 
most recent MDI observations processed are from early 2009. Overlaid are contours of the gross 
magnetic longitudinal magnetic field strength. The branch of the torsional oscillation pattern corre- 
sponding to Cycle 24 is well established and clearly visible, but is moving more slowly towards the 
equator than did the corresponding branch in Cycle 23. The low-latitude near-surface flow profile 
in this data set is best matched by that at date 1997.3, as shown by the left-hand vertical line in the 
figure. At this point in cycle 23, substantial magnetic activity was just starting, as illustrated by the 
latitudinal magnetic profiles in Figure [2] The level of activity seen for the current cycle is much 
lower, but what hints of new-cycle activity there are seem to be at the correct latitude. Note also 
that there seem to have still been a few old-cycle active regions around at this point in early 1997; at 
present only a few very weak old-cycle regions are seen, for example, NOAA 1 1016 on 2009 April 
29 — May 1. To give an idea of the depth dependence of the flows, Figure Q] also shows twelve- 
month averages of the residuals in the r, 9 plane for four selected low-activity epochs. The radial 
detail of these profiles should not be over-interpreted, due to the poorer signal-to-noise ratio of the 
inversions in the interior and the varying resolution of different inversions, but the penetration of 
the flow belt into the bulk of the convection zone is clear. 

How does the evolution of the flows during the current minimum correspond to what was 
seen in the previous one? To quantify this, we consider the linear correlation between the near- 
surface residuals 80(9, r,t±) at each time step t\ and those at a reference epoch t , 80,(9, r,t ). 



Available from |http : / /nsokp . nso . edu/dataarch . html| 
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The correlation coefficient between two variables x(t ) and x(t\) is defined (lBevingtonll 19691) as 



1 °' lj {£J^o) - x(t )] 2 } 1/2 {EJ^(ii) - ^i)] 2 } 1/2 ' 1 ; 

where in this case = 5Cl(Ti,t) and the index i enumerates inversion mesh points with radius 
0.97 < r / R Q < 0.995 and < 9 < 45deg. In Figure [3] we show the results of this analysis for 
four selected reference epochs. For a reference epoch of 1995.5 (Figure^), about a year before 
the last solar minimum, the correlation reaches a maximum value around the second half of 2006 
to early 2007 and then declines sharply. When the reference epoch is at 1996.5, corresponding 
to the previous solar minimum, (Figured) the correlation coefficient reaches a maximum around 
the middle of 2008, and then starts to decline. A reference epoch of 2006.5 (Figure |3k) gives a 
maximum correlation at the beginning of the observations, 1995.5, while a reference epoch of the 
most recent set of observations, 2009.2, gives a correlation peaking at 1997.3, about eleven months 
after solar minimum (Figure [3]i). 

To study the temporal evolution of the apparent cycle length, for each data set from the be- 
ginning of 2007 onwards we evaluate equation (QQ) over latitudinal mesh points up to 45 deg, but 
separately at each radial mesh point from 0.973 to 0.995-R©, and find the date of the highest cor- 
relation value in the early part of the cycle. The difference between this date and the reference 
epoch gives an estimate of cycle length in years; we then take the mean and standard deviation 
of these lengths at each time step, thus obtaining an estimate of the effective cycle length and its 
uncertainty as a function of time. The results, seen in Figure |4l seem to suggest that the effective 
length of the cycle measured in this manner has been slowly increasing through 2007 - 2008. 

An alterna t ive approach to describing the temporal variation of the flows, in troduced by 



Voro ntsov et al.l (|2002|) . is to fit the variation at each location with a sinusoid function. iHowe et al 
(|2005|) showed that, for the data through 2004, the variation was well described by fitting with an 
11 -year period sine wave and its second harmonic, but did not find it practical to investigate other 
possible cycle lengths with only nine years of observations. With nearly fourteen years of data, we 
can now test the fit with other cycle lengths. The time variation of the rotation at each location is 
modeled as 

fifr, 0, t) = tt (r, 6) + Ai sin(27rf/P) + A 2 cos(2nt/P) + A 3 sm(Ant/P) + A 4 cos(4vrt/P) , (2) 

where P is the period in years and time t is also measured in years. We then vary P over the range 
from 9 to 14 years and find the value which gives the lowest x 2 at each mesh point, before taking 
the mean result and its standard deviation over the range 0.97 < r/R Q < 0.995. Again, we see 
a tendency to gradually increasing perceived cycle length in these best-fit periods, as shown by 
the curve in Figure HI Taking these two analyses together, we estimate that the length of the cycle 
in the region equatorwards of 45 deg increased from about 1 1.2 to 12.2 years between early 2007 
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and early 2009. We have carried out a similar sine-wave analysis at higher latitudes; poleward of 
60 deg the cycle length remains almost constant at about 12 years over the same period, while at 
intermediate latitudes the behavior is more complicated and not clearly periodic. In other words, 
the flow configuration, in agreement with other indicators, seems to indicate a prolonged minimum, 
with the flow belts (and, presumably, the almost-undetectable magnetic activity belts) moving more 
slowly towards the equator than was the case during the previous minimum. 



4. CONCLUSION 

The zonal flow pattern during the current solar minimum, considered alongside that for the 
previous minimum and the magnetic butterfly diagram, suggests that the minimum was reached 
during 2008, giving a length of approximately 12 years for Cycle 23. The flow band associated 
with the new cycle has been moving more slowly towards the equator than that observed during 
the previous minimum. 

It would not be surprising to see a rapid increase in activity in the near future; it is tempting 
to suggest that if the current low levels persist for much longer we may indeed be looking at an 
unusually weak cycle. 

In follow-up work, we will investigate the behavior of the flows in more detail, including the 
poleward branch and the flow in the deeper layers. 

This work utilizes data obtained by the Global Oscillation Network Group (GONG) pro- 
gram, managed by the National Solar Observatory, which is operated by AURA, Inc. under a 
cooperative agreement with the National Science Foundation. The data were acquired by instru- 
ments operated by the Big Bear Solar Observatory, High Altitude Observatory, Learmonth Solar 
Observatory, Udaipur Solar Observatory, Instituto de Astrofisica de Canarias, and Cerro Tololo In- 
teramerican Observatory. The Solar Oscillations Investigation (SOI) involving MDI is supported 
by NASA grant NAG 5-13261 to Stanford University. SOHO is a mission of international co- 
operation between ESA and NASA. RK, and RH in part, were supported by NASA contracts S- 
92698-F, NAG 5-1 1703, and NNG 05HL41I. NSO/Kitt Peak data used here were produced coop- 
eratively by NSF/NOAO, NASA/GSFC, and NOAA/SEL; SOLIS data are produced cooperatively 
by NSF/NSO and NASA/LWS. 
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Fig. 1. — (a) Rotation-rate residuals at O.99i? from MDI and GONG. Overlaid contours show 
the gross longitudinal magnetic field strength from KPVT/SOLIS, at 5G intervals. The leftmost 
solid vertical light-blue line shows the date, 1997.3, at which the low-latitude flow configuration 
best matches that in the most recent (2009.2) data set, and rightmost vertical line the date, 2006.4, 
where it best matches that in the earliest data set (1996.5), while the horizontal lines show the 
respective location of the flow bands and the slanted lines schematically indicate the migration of 
the equatorward branch. The lower panels show 12-month averages of the rotation-rate residuals 
in the r, plane for epochs starting at (b) 1995.5, (c) 1996.3, (d) 2006.5, (e) 2008.2. 
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Fig. 2. — KPVT/SOLIS gross longitudinal magnetic field strength as a function of latitude, for the 
2009.2 epoch of the most recent zonal flow observations (solid line) and for the epoch, 1997.3, 
during the previous minimum where the flows most closely match the most recent ones (dashed 
line). 
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Fig. 3. — Correlation between near-surface rotation residuals at each time sample and those at 
selected reference epochs. Triangles represent GONG data and circles MDI. The reference epochs 
are: (a) one of the earliest sets of GONG observations at 1995.5; (b) the solar minimum at 1996.5; 
(c) the epoch 2006.5, chosen to give the best match to the earliest observations; and (d) the last set 
of observations we have at 2009.2. 
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Fig. 4. — Estimated length of cycle as a function of last observing time, averaged over 0.97 < 
r/ Rq < 0.995 and latitude less than 45 deg. Results are derived from time lag for best correlation 
(symbols, with the error bar representing the standard deviation of the mean), and length of cycle 
estimated by sinusoid fit (solid line, with standard deviation of the mean shown by the dotted lines). 



